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We analyze the indirect exchange interaction between two two-state systems, e.g., spins 1/2, 
subject to a common finite-temperature environment modeled by bosonic modes. The environmental 
modes, e.g., phonons or cavity photons, are also a source of quantum noise. We analyze the coherent 
vs noise-induced features of the two-spin dynamics and predict that for low enough temperatures the 
induced interaction is coherent over time scales sufficient to create entanglement. A nonperturbative 
approach is utilized to obtain an exact solution for the onset of the induced interaction, whereas 
for large times, a Markovian scheme is used. We identify the time scales for which the spins 
develop entanglement for various spatial separations. For large enough times, the initially created 
entanglement is erased by quantum noise. Estimates for the interaction and the level of quantum 
noise for localized impurity electron spins in Si-Ge type semiconductors are given. 
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I. INTRODUCTION 

The idea that exchange of fermionic or bosonic excita- 
tions can lead to physical interactions in solid state is not 
new.^ Recently, such induced interactions have received 
attention due to the possibility to experimentally observe 
quantum dynamics in nanoscale devices i^i^i^i^i^i^ii^ii^ In 
this work we explore the dynamics of two qubits (two 
two-state quantum systems), e.g., electron spins 1/2, 
placed a distance d apart, as they are entangled by 
common thermalized bosonic environment without direct 
spatial electron wave function overlap. At the same time 
quantum noise originating from the same bosonic field 
(e.g., phonons) ultimately erases the generated entan- 
glement for large enought times. We demonstrate that 
the indirect exchange interaction induced by the bosonic 
thermal fieldi^ can, in some cases, be comparable with 
other qubit-qubit couplings. 

Extensive studies have been reportediiii^ii^ii^ of 
the decay of quantum correlations between qubits 
subject to individual (local) or common environ- 
mental noise. In the presence of quantum noise, 
entanglement was shown to decay very fast and, 
in some cases, vanish at finite timesJ^ii^ On the 
other hand, the idea that common bosonic as well 
as fermionic environment is able to entangle the 
qubits has also been advanced ^iiO^ii^aiis 09,20^^^ 

For fermionic environment, this effect has been at- 
tributed to Rudermann-Kittel-Kasuya-Yosida (RKKY) 
type interaction a-'^^i-'^^'^^i^^i^^i^^ 'Si? and it has recently been 
experimentally demonstrated for two-qubit semiconduc- 
tor nanostructures^i^ 

In this work, we investigate the two competing phys- 
ical effects of a thermalized bosonic environment (bath) 



in which two qubits are immersed. Specifically, the in- 
duced interaction, which is effectively a zero-temperature 
effect, and the quantum noise, originating from the same 
bath modes, are derived within a uniform treatment. 
We study the dependence of the induced coherent vs 
noise (decoherence) effects on the parameters of the bath 
modes, the qubit system, and their coupling, as well as 
on the geometry. Specific applications are given for spins 
interacting with phonons in semiconductors. 

We assume that the spins are identically coupled with 
the modes of a thermalized bosonic bath, described by 
Hb = X^k J "^k.^Ok fOk.^ , where we introduced the polar- 
ization index ^ and set h = 1. With qubits represented 
as spins 1/2, the external magnetic field is introduced in 
the qubit Hamiltonian 



?)/2, 



(1.1) 



as the energy gap A between the up and down states 
for spins 1 and 2, with the spins labeled by the super- 
scripts. A natural example of such a system are spins of 
two localized electrons interacting via lattice vibrations 
(phonons) by means of the spin-orbit interaction>i i^^i^^i^^ 
Another example is provided by atoms or ions in a cav- 
ity, used as two-state systems interacting with photons. 
For each type of phonon/photon, the interaction will be 
assumediiS£iS2i2^ linear in the bosonic variables. The 
spin-boson coupling for two spins j ~ 1,2 has the form 



H, 



SB 



E E ' 

J — 1,2 m—x^y.z 



where a-L are the standard Pauli matrices and 



XL 



e 



(ak,c + alk.e) 



(1.2) 



(1.3) 
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Here, as before, the index ^ accounts for the polarization 
and Tj is the position of jth spin. The overall system 
is described by the Hamiltonian H = Hs + Hb + Hsb- 
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Our emphasis will be on comparing the relative impor- 
tance of the coherent vs. noise effects of a given bosonic 
bath in the two-qubit dynamics. We do not include other 
possible two-qubit interactions in such comparative cal- 
culation of dynamical quantities. 

For the analysis of the induced exchange interaction 
and quantum noise for most time scales of relevance, 
it is appropriate to use the Markovian approach22i2^i2£ 
that assumes instantaneous rethermalization of the bath 
modes. In Section [ITl we derive in a unified formulation 
the bath-induced spin-spin interaction and noise terms 
in the dynamical equation for the spin density matrix. 
The onset of the interaction for very short times is also 
investigated within an exactly solvable model presented 
in Sections IIIII and HVl Specifically, in Section HVl we dis- 
cuss the onset and development of the interaction Hamil- 
tonian, as well as the density matrix. It is shown that 
the initially unentangled spins can develop entanglement. 
We find that the degree of the entanglement and the time 
scale of its ultimate erasure due to noise can be controlled 
by varying several parameters, as further discussed for 
various bath types in Section |Vl Estimates for the in- 
duced coherent interaction in Si-Ge type semiconductors 
are presented in Section IVll The coherent interaction in- 
duced by phonons is, expectedly, quite weak. However, 
we find that in strong magnetic fields it can become com- 
parable with the dipole-dipole coupling. 



II. COHERENT INTERACTION AND 
QUANTUM NOISE INDUCED BY 
THERMALIZED BOSONIC FIELD 

In this section we present the expressions for the in- 
duced interaction and also for the noise effects due to 
the bosonic environment, calculated perturbatively to the 
second order in the spin-boson interaction, and with the 
assumption that the environment is constantly reset to 
thermal. Specific applications and examples are consid- 
ered towards the end of this section, as well as in Sec- 
tions |V] and Ell 

The dynamics of the system can be described by the 
equation for the density matrix 



ip{t) = [H,p{t)]. 



(2.1) 



In order to trace over the bath variables, we carry 
out the second-order perturbative expansion. This dy- 
namical description is supplemented by the Markovain 



assumption- 



27.28.29.30 



of resetting the bath to thermal 



equilibrium, at temperature T, after each infinitesimal 
time step, as well as at time t = Q, thereby decoupling 
the qubit system from the environment 3^1^ This is a 
physical assumption appropriate for all but the shortest 
time scales of the system dynamics ^ii^i^ It can also 
be viewed as a means to phenomenologically account in 
part for the randomization of the bath modes due to their 
interactions with each other (anharmonicity) in real sys- 
tems. This leads to the master equation for the reduced 



density matrix of the qubits ps{t) = Trspit), 

ips{t) = [Hs,ps{t)] (2.2) 

OO 

- I j dt'TrB[HsB,[HsB{t' -t),pBPs{t)]], 



where i?ss(r) = e^(^B+^s)^i/sBe-^(^«+^s)T 



PB = 



•^ala^/knT^ and Z = 1/11^(1 - e-'^"/''^^) is 



the partition function. Analyzing the structure of the 
integrand in Eq. (|2.2p . one can obtain the equation with 
explicitly separated coherent and noise contributions, see 
Appendix \^ 



ipsit) = [H,s,ps{t)]+iMps{t). 
Here the effective coherent Hamiltonian H^s is 



(2.3) 



m—x.y,z 

(2.4) 

The expressions for the amplitudes x™(d), x™(d), 
77™ (d), and 77™ (d) will be given shortly. The first three 
terms following Hg constitute the interaction between 
the two spins. We will argue below that the leading in- 
duced exchange interaction is given by the first added 
term, proportional to x™(d). The last term corrects the 
energy gap for each qubit, representing their Lamb shifts. 

The explicit expression for the noise term is very cum- 
bersome. It can be represented concisely by introducing 
the noise superoperator Af , which involves singlc-qubit 
contributions, which are usually dominant, as well as 
two-qubits terms 



M 



(2.5) 



where the summations are over the components, m = 
X, y, z, and the qubits, z, j = 1, 2. The quantities entering 
Eq. (|2.5p can be written in terms of the amplitudes x™(d), 
77™ (d), and 77™ (d), in a compact form, by introducing the 
superoperators La{0i)02 = {Oi,02}, ^(Oi, 02)03 = 
O1O3O2, and L±(Oi, O2) = L(Oi, O2) ± L(02, Oi), 



M^"(d) = 77r(d)[2L«„,a^)-L,(a;,0 
+ '7r(d) [L+(aj„,0-i„(aj„?/, 

- 7xr(d)i-«„^4), 
where we defined q^^ — ^[al, cr^J and 

Mm = 77r(0)[2L«,O-L,(a>^ 

+ 77r(o)i+«„^^) 

- *xr(0)U-(a^,C)+ia«?;( 



(2.6) 



(2.7) 
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The amplitudes in Eqs.^M EJ]), calculated for 
the interaction defined in Eqs. ljf .21 II. 3p . are 



Vdk 



,2 [Jk,4 COS (k • d) 



(2.8) 
(2.9) 



^ — oo 

X ^ (5 (tJk,4 + (1 - '5™.2)gA), 

g=±l 



and 



xTid) = -(l-<5™,.)^ 7Z^|g™^|2cos(k.d) 
^ J (27r) 

(2.10) 



q=±l ^ 



r;r(d) = (1-5™,.) 



E 



(2^ 



' coth 



(2.11) 
wk,^ Acos(k-d) 



2k^T Lo- 



A2 



Here the principal values of integrals are assumed. 

Note that x™(d) appears only in the induced inter- 
action Hamiltonian in Eq. (|2.4p . whereas 7y™(d), xr(d), 
and f7™(d) enter both the interaction and noise terms. 
Therefore, in order to establish that the induced inter- 
action can be significant for some time scales, we have 
to demonstrate that (d) can have a much larger mag- 
nitude than the maximum of the magnitudes of r7™(d), 
X™(d), and rj^{d). The third and fourth terms in ex- 
pression for the interaction (|2.4p are comparable to the 
noise and therefore have no significant contribution to 
the coherent dynamics. 

Because of the complexity of the expressions for the 
noise terms within the present Markovian treatment, in 
this section we will only compare the magnitudes of the 
coherent vs noise effects. In the next section, we will 
discuss a different model for the noise which will allow a 
more explicit investigation of the time dependence. 

For the rest of this section, we will consider an illus- 
trative one-dimensional (ID) example favored by recent 
experimentsj^i^ leaving the derivations for higher dimen- 
sions to Sections |V] and IVII We comment that the ID 
geometry is also natural for certain ion-trap quantum- 
computing schemes, in which ions in a chain are subject 
to Coulomb interaction, developing a variety of phonon- 
mode lattice vibrationsi ^^i'^^'^^ 

In ID geometry we allow the phonons to propagate 
in a single direction, along d, so that k • d = k\d\. 
Here, for definiteness, we also assume the linear disper- 
sion, ojfe — Cgk, since the details of the dispersion re- 
lation for larger frequencies usually have little effect on 



the decoherence properties. Furthermore, we ignore the 
polarization, gj^^ g™'{uj). Another reason to focus 
on the low-frequency modes is that an additional cutoff 
LUc resulting from the localization of the electron wave 
functions, typically much smaller than the Debye fre- 
quency, will be present due to the factors g'^"-{uo). The 
induced interaction and noise terms depend on the am- 
plitudes xr(d), ?7"'(d), xr(d), and ?7r(d), two of which 
can be evaluated explicitly for the ID case, because of 
the 5-functions in Eas. (|2.91 I2.10p . However, to derive 
an exphcit expression for xr(d) and ■q^{d), one needs 
to specify the w-dependence in \g™-{u!)\^. For the sake 
of simplicity, in this section we approximate [(/"^(w)!^ by 
a linear function with superimposed exponential cutoff. 
For a constant ID density of states T(w) = V/2ttCs, this 
is the Ohmic-dissipation condition, i.e.. 



^uj"" exp(— w/wc). 



(2.12) 



with n = 1 (the case when n > 1 is considered in Section 

In most practical applications, we expect that A ^ Wc- 
With this assumption, we obtain 



1 + {UJ, \d\/CsY 



xr(d) = 



c- X A Aldl 

(1 - Ojrt,z) COS 

tOc Cs 



(2.13) 



(2.14) 



and 



(d) = aTu;,^ (1 - <5„,,) - coth 

2 Ll)c 2kb-L 



Aldl 



cos ■ 



(2.15) 

The expression for »7s(d) could not be obtained in closed 
form. However, numerical estimates suggest that ?7s(d) is 
comparable to ?7™(d). At short spin separations d, ?7s(d) 

is approximately bounded by — q;™A In ^ exp(— ^r^), 
while at lager distances it may be approximated by 
a™A^ sin coth t • "^^^ level of noise may be es- 
timated by considering the quantity 



M=maxhr(d),xr(d),77r(d) 

|d| 



The interaction Hamiltonian takes the form 



(2.16) 



This induced interaction is temperature independent. It 
is long-range and decays as a powerlaw for large d. For 
the super-Ohmic case, n > 1, one obtains similar be- 
havior, except that the interaction decays as a higher 
negative power of d, as will be shown in Section IVl 

If the noise term were not present, the spin system 
would be governed by the Hamiltonian Hg + Hi^t- To 
be specific, let us analyze the spectrum, for instance. 
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FIG. 1: Short-time correction to the induced exchange inter- 
action for the Ohmic case. The arrow shows the order of the 
curves for increasing [iJc|d|/cs = 0, 1, 10. 



for 



a\, af = 0. The two-qubit states consist of 



the singlet (|ti) - Ut>)/%/2 and the spHt triplet |tT), 
(ITi) + UT»/V2, and lii), with energies E2 = -4Xc, 
Eo = —A, El — 4xc, E3 = A, respectively. The en- 
ergy gap \Ei — £2] between the two entangled states is 
defined by 4Xc (= 4Xc)- In the presence of noise, the os- 
cillatory, approximately coherent evolution of the spins 
can be observed over several oscillation cycles provided 
that 2aY'uJc/[l + (wcldl/c,)^] > M. The energy levels 
will acquire effective width due to quantum noise, of or- 
der ?7™(d). This interplay between the interaction and 
noise effects is further explored within an exactly solvable 
model in the next section. 



III. EXACTLY SOLVABLE MODEL 

In this section, we consider a model appropriate for 
short times which does not invoke the Markovian 
assumption of the rethermalization of the bath modes. 
This model is particularly suitable for investigating the 
onset of the system's dynamics. While the noise effects 
are quantitatively different in this model, the qualitative 
interplay of the coherent and noise effects in the dynamics 
is the same as in the Markovian approach. Furthermore, 
we will show that the induced interaction is consistent 
with the one obtained within Markovian approach in the 
previous section. 

We point out that, due to the instantaneous rethermal- 
ization assumption (resetting the density matrix of the 
bath to thermal), in the Markovain formulation it was 
quite natural to assume that the bath density matrix is 
also thermal at time t = 0; the total density matrix re- 
tained an uncorrelated-product form at all times. In the 
context of studying the short-time dynamics, in this sec- 
tion the choice of the initial condition must be addressed 
more carefully. In quantum computing applications, the 
initially factorized initial condition has been widely used 
for the qubit-bath systeiri^i^i^ 
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FIG. 2: Development of the concurrence as a function of time, 
calculated with = 1/20 and ksT/ujc = 1/20. The curves 
correspond to various spin-spin separations, as can be read off 
the left inset, which shows the distribution of the concurrence 
in the \d\-t plane. The right inset presents the dynamics of 
the diagonal density matrix elements P|| = (IT \ps\ TT)> etc., 
on the same time scale. 



This choice allows comparison with the Markovian re- 
sults, and is usually needed in order to make the short- 
time approximation schemes tractable^^ specifically, it 
is necessary for exact solvability of the model considered 
in this and the next sections. 

A somewhat more "physical" excuse for factorized ini- 
tial conditions has been the following. Quantum compu- 
tation is carried out over a sequence of time intervals dur- 
ing which various operations are performed on individual 
qubits and on pairs of qubits. These operations include 
control gates, measurement, and error correction. It is 
usually assumed that these "control" functions, involv- 
ing rather strong interactions with external objects, as 
compared to interactions with sources of quantum noise, 
erase the fragile entanglement with the bath modes that 
qubits can develop before those time intervals when they 
are "left alone" to evolve under their internal (and bath 
induced) interactions. Thus, for evaluating relative im- 
portance of the quantum noise effects on the internal (and 
bath induced) qubit dynamics, which is our goal here, we 
can assume that the state of qubit-bath system is "reset" 
to uncorrelated at t ^ 0. 

It turns out that the resulting model is exactly solvable 
for the Zeeman splitting A = 0, and provided that only 
a single system operator enters the expression (|1.2p for 
the interaction. Here we take = = 0, while ^ 
0. We derive the exact solution and demonstrate the 
emergence of the interaction (|2.17p . 

With the above assumptions, one can utilize the 
bosonic operator techniques'^ to obtain the reduced den- 
sity matrix for the system (|1.21 II. 3p , 



A,A' 



(3.2) 



piO) = Ps{0)pb{0). 



(3.1) 



where the projection operator is defined as Pa = 
IA1A2) (A1A2I, and |Aj) are the eigenvectors of cr^. The 



5 



exponent in Eq. (|3.2p consists of the real part, which leads 
to decay of off-diagonal density-matrix elements resulting 
in decoherence, 



Re2xx'{t) - -J2Gk{t,T) \{X[ - X^f + (A^, - A2) 



2 cos 



/ ujk |d| 



)(A'i-Ai)(A'2-A2) 



(3.3) 



and the imaginary part, which describes the coherent 
evolution. 



Im£AA'(t) -Jl^fcW 



(A1A2 — X[X'2)- 



Here we defined the standard spectral functions^i^ 



(3.4) 



G,(t,T) = 2^sin^(^^Jcoth 



2kBT 



(3.5) 



F{t) - a>;^/(wcO"- In Figured we plot J^{t), defined 
via F{t) — J^{t)al.a'^, for the Ohmic case as a function of 
time for various spin-spin separations. 

Let us now explore the role of the decoherence term 
(13. 3p . In the exact solution of the short-time model, 
the bath is thermalized only initially, while in the per- 
turbative Markovian approximation, one assumes that 
the bath is reset to thermal after each infinitesimal time 
step. Nevertheless, the effect of the noise is expected 
to be qualitatively similar. Since the short-time model 
offers an exact solution, we will use it to compare the 
coherent vs noise effects in the two-spin dynamics. We 
evaluate the concurrence j^^*^ which measures the entan- 
glement of the spin system and is monotonically related 
to the entanglement of formation. '^^■'"^ For a mixed state 
of two qubits, ps, we first define the spin-flipped state, 
Ps — ^y^y P*s ^y^y^ ^^d then the Hermitian operator 

= \/ y/psPS\/ps, with eigenvalues Ai=i,2,3,4. The con- 
currence is then given™ by 



and 



Ck{t) 



\9k\ 



(ujkt - sinwfct) . 



(3.6) 



Calculating the sums by converting them to integrals 
over the bath- mode frequencies uj in Eqs. p.3p and (|3.4p . 
assuming the Ohmic bath n — I for T > one ob- 
tains a linear in time t large-time behavior for both the 
temperature-dependent real part and for the imaginary 
part. The coefficient for the former is ^kT, whereas for 
the latter it is ~ tOc- For super-Ohmic models, n > 1, 
the real part grows slower, as was also noted in the 
literature,22i^i2^ 

Let us first analyze the effect that the imaginary part 
of £aa' (t) has on the evolution of the reduced density 
matrix, since this contribution leads to the induced in- 
teraction. If the real part were not present, i.e., omitting 
Eq.(IS31) from Eqs.dST!]), dSH), and dSl]), we would ob- 
tain the evolution operator in the form exjp[—iHintt — 
iF{t)t]. The interaction iJi„t comes from the first term 
in Eq. 



Hint — ' 



2<r(n)c>,^ 



n/2 



n arctan 



(3.7) 



This expression is the same as the results obtained within 
the Markovian scheme, cf., Eas. (|TT7|) . (j^ . and Sec- 
tion |Vl The operator F{t) is given by 



2 

X^ X 



, D(uj) \q(iij)\ sinujt 
duj ^ cos 

UJ LUt 



(3.8) 

It commutes with H-mt and therefore could be viewed as 
the initial time-dependent correction to the interaction. 
In fact, it describes the onset of the induced coherent 
interaction; note that F{0) = —Hint, but for large times 



Cips) 



max 



0, 2 max Xi — Xj 



(3.9) 



In Figure [2l we plot the concurrence as a function of 
time and the spin-spin separation, for the (initially un- 
entangled) state |tT): and n = 1. One observes decaying 
periodic oscillations of entanglement. We should point 
out that the measure of entanglement we use here — the 
concurrence — provides the estimate of how much entan- 
glement can be constructed provided the worst possible 
scenario for decomposing the density matrix is realized; 
see Refs. 113, HSi for details and definitions. Therefore, 
one expects the entanglement that one can make use of 
in quantum computing to be no smaller than the one 
presented in Figure [2] In the next section, additional 
quantities are considered, namely, the density matrix el- 
ements, which characterize the degree of coherence in the 
system's dynamics. 



IV. ONSET OF THE INTERACTION AND 
DYNAMICS OF THE DENSITY MATRIX 

Let us now investigate in greater detail the onset of 
the induced interaction, the time-dependence of which 
is given by F{t). In Figure [1] we have shown the magni- 
tude of F{t), as a function of time for various qubit-qubit 
separations and n = 1. The correction is initially non- 
monotonic, but decreases for larger times as mentioned 
above. The behavior for other non-Ohmic regimes is ini- 
tially more complicated, however the large time behavior 
is similar. 

It may be instructive to consider the time dependent 
correction Hpit) to the interaction Hamiltonian dur- 
ing the initial evolution, corresponding to F{t). Since 
Fit) commutes with itself at different times, as well as 
with -ffjnt: it generates unitary evolution according to 
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FIG. 3: The magnitude of the time-dependent Hamiltonian 
corresponding to the initial correction as a function of time 
and distance. The Ohmic (n = 1) case is shown. The inset 
demonstrates the onset of the cross-qubit interaction on the 
same time scale. 



exp[ 



-i J* dt'Epit')], with Hpit) = d[tF{t)]/dt, 



therefore 



Hpit) = a>2«„r(n) (4.1) 

X [u{uJc\d\/Cs - LOct) + u{uJc\d\/Cs + UJct)], 

where u(^) = cos[n arctan(C)]/[l + ^^1"/^ The above 
expression is a superposition of two waves propagating in 
opposite directions. In the Ohmic case, n ~ 1, the shape 
of the wave is simply u{^) = 1/(1 + f^). In Figure [31 
we present the amplitude of Hp{t), defined via Hp{t) = 
Hp'^l'^x' ^ "^^ll ^ the sum of ifint and Hp{t), for n = 1. 
One can observe that the "onset wave" of considerable 
amplitude and of shape m(^) propagates once between 
the qubits, "switching on" the interaction. It does not 
affect the qubits once the interaction has set in. 

To understand the dynamics of the qubit system and 
its entanglement, let us again begin with the analysis 
of the coherent part in Eq. p.2p . After the interaction, 
i/int, has set in, it will split the system energies into two 
degenerate pairs Eg = Ei = —Hint and E2 — E^ — Hint- 
The wave function is then \tp{t)) = exp(— ii/intOI'0(O))- 
For the initial "up-up" state, \ip{0)) — |TT), it develops 
as IV'(i)) = ITT)cos7iintt + isinHintt, where Hi^t = 
Hint'^m'^m- O^i^ Can easily notice that at times tp — 
7r/4Hint, 37r/4Hint, ■ • maximally entangled (Bell) states 
are obtained, while at times to ~ 0, 7r/27iint, "'/'^inti ■ • 
the entanglement vanishes; these special times can also 
be seen in Figure O The coherent dynamics obtained 
with the Markovian assumption is the same. 

However, the coherent dynamics just described is only 
approximate, because the bath also induces decoherence 
that enters via Eq. p.3p . The result for the entanglement 
is that the decaying envelope function is superimposed 
on the coherent ocsillations described above. The mag- 
nitudes of the first and subsequent peaks of the concur- 
rence are determined only by this function. As temper- 



ature increases, the envelope decays faster resulting in 
lower values of the concurrence; see the inset in Figure [H 
Although in the Markovian approach, presented in Sec- 
tion [Til the noise is quantitatively different, one expects 
qualitatively similar results for the dynamics of entangle- 
ment. 

Note also that nonmonotonic behavior of the entan- 
glement is possible only provided that the initial state 
is a superposition of the eigenvectors of the induced in- 
teraction with more than one eigenvalue (for pure initial 
states); see Eqs. p.2p - (l3.4p . For example, taking the ini- 
tial state (|iT) + |n))/\/2, in our case would only lead to 
the destruction of entanglement, i.e., to a monotonically 
decreasing concurrence, similar to the results of Refs. [ij, 

m 

For the model that allows the exact solution, i.e., for 
Hs ~ 0, one can notice that there is no relaxation by 
energy transfer between the system and bath. The ex- 
ponentials in Eas. (|3.2p . with (|3.3p . suppress only the off- 
diagonal matrix elements, i.e., those with A 7^ A'. It hap- 
pens, however, that at large times the d-dependence is 
not important in Eq. (j3.3p . and Ile£,xy{t — > 00) vanishes 
for certain values of A 7^ A'. In the basis of the qubit-bath 
interaction, cr^cr^, the limiting t — s- 00 density matrix for 



our initial state (|tT)) is 
which takes the form 



pit 



/ 3 00-1 
0110 
0110 
V -1 3 



(4.2) 



in the basis of states |tT), Iti); lit), and The sig- 

nificance of such results, see also Ref. lUi, is that in the 
model with Hs = and nonrethermalizing bath not all 
the off-diagonal matrix elements are suppressed by deco- 
herence, even though the concurrence of this mixed state 
is zero. 

The probabilities for the spins to occupy the states 
ITT), ITi), liT), and III) are presented in Figure[31 For the 
initial state |TT)i only the diagonal and inverse-diagonal 
matrix elements are affected, and the system oscillates 
between the two states |TT) and |||), as mentioned earlier 
in the description of the coherent dynamics, while deco- 
herence dampens these oscillations. In addition, decoher- 
ence actually raises the other two diagonal elements to a 
certain level, see Ea. (l4.2p . The dynamics of the inverse- 
diagonal density matrix elements is shown in Figure [51 



V. OHMIC AND SUPER-OHMIC BATH 
MODELS IN GENERAL DIMENSION 

Let us generalize the results of the previous sections ob- 
tained primarily for the Ohmic bath model and ID geom- 
etry. In Section[TTl we considered the ID case with Ohmic 
dissipation with the Markovian approach. In the general 
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FIG. 4: Dynamics of the occupation probabilities for the 
states lit), III), lit), and |J,|,). The parameters are the same 
as in Figure [2] The inset shows the structure of the reduced 
density matrix (the nonshaded entries are all zero). 
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the same system as in Figure ID Note that Impn^n = 0. 



FIG. 6: The magnitude of the induced exchange interaction 
Hamiltonian for the Ohmic rim ~ 1 and super-Ohmic rim = 
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FIG. 7: The magnitude of the noise for the Ohmic rim = 1 
and super-Ohmic n,„ = 2, 3, 4, 5 baths, for 2k-BT /uoc ~ 0.01. 



case, let us consider the Markovain model and, again, as- 
sume that A/ujc is small. We will also assume that the 
absolute square of the mth component of the spin-boson 
coupling, when multiplied by the density of states, can 
be modelled by exp(-a;/a;c); see Eq. ([21121). The 

integration in Eq. (|2.8p can then be carried out in closed 
form for any Um — 1,2,... . The induced interaction 
(|2.17p is thus generahzed to 

i?int = - E <„^c'"'^™^™ (5.1) 
m—x,y,z 



2T{n 


™)Re(l + iu!c\d 






l + {uj,\d\/c,f 





With the appropriate choice of parameters, the result 
for the induced interaction, but not for the noise, coin- 
cides with the expression for iJint obtained within the 
short-time model. From Eq. (|5.ip one can infer that the 
effective interaction has the large-distance asymptotic be- 
havior |d| for even rt,„, and |d| "™ ^, for odd rim- 

Let us now consider the noise terms. The amplitudes 
(|2.91 12.10p . entering the decoherence superoperator, de- 



pend on Um only via the prefactor (A/cjc)"" ■ The re- 
maining amplitude (|2.1ip has to be calculated numeri- 
cally. The spin-spin separation up to which the "coher- 
ent" effects of the induced interaction can be observed in 
the time-dependent dynamics, is defined by comparing 
the magnitudes of the interaction and noise amplitudes. 
It transpires that the comparison of the amplitude in 
Ea. (j2.17p vs Eq. (|2.16p for rim — 1, also suffices for an 
approximate estimate for the super-Ohmic long 
as 77™ (d) does not dominate in M. In Figure [HI we plot 
the amplitudes of Ea. (|5.ip . i.e., Hint = Tr{Hint(Tm<Tm) / ^1 
for different values of n^. Figure [7] shows the magni- 
tude of the noise super-operator. Figures 6 and 7 sug- 
gest that the noise amplitudes can be made sufficiently 
small with respect to the induced interaction, at small 
values of A/wc, for any n^. We note that for quantum 
computing applications one usually assumes the regime 
ksT <^ A <^ LOc- The temperature dependence in Fig. 7 
becomes insignificant for fc^T/A <^ 1, which is approx- 
imately to the right of A/cjc ^ 0.01. To the left of 
A/ujc ~ 0.01, by reducing the temperature one can fur- 
ther reduce the values of the noise amplitudes even for 
small A. 

In higher dimensions the structure of 



8 



in the k-space becomes 
tJk.^ is nearly isotropic, 



will 



important. Provided 
the integrals entering 
include (in 3D) a factor 

2 

which can be 



27r 7T ^ 

/o '^'^ lo si'^ ^ 5^9.^,^ COS (fc |d| COS ( 
written as [fl''{u;,k\d\) - /2™(w, fc|d|)9/a|d|] cos(/c|d|), 
e.g., Eqs. (|B.21 IB.3P in Appendix |B] When the depen- 
dence of Z™;/™ on A:|d| is negligible, the interaction 

is simply Hint i/int|{„„}^a - id/d\d\)Hint\{n^}->b, 

where a and b are sets of three integers representing 
the a;-dependence of /™(cj,fc|d|) and f^{uj,k\d\). 
Otherwise, a more complicated dependence on |d| is 
expected. The noise superoperator can be treated 
similarly. As a result one can see that the form of 
the interaction depends more on the structure of g^^^ 
than on the dimensionality via the phonon density of 
states. However, reduced dimensionality in the k-space 
might allow for better control over the magnitude of 
the interaction by external potentials that modify the 
spin-orbit coupling. 



VI. PHONON INDUCED COHERENT 
SPIN-SPIN INTERACTION VS NOISE FOR 
P-DONOR ELECTRONS IN Si AND Ge 

As a specific example, let us consider a model of two 
localized impurity-electron spins of phosphorus donors in 
a Si-Ge type semiconductor, coupled to acoustic phonon 
modes by spin-orbit interaction. In what follows, we 
first obtain the coupling constants g^^^, entering Eq. (|1.3p . 
which define the interaction and noise amplitudes. A 
brief discussion is then offered on the possibility of hav- 
ing an Ohmic bath model realized in ID channels. The 
rest of the present section is devoted to calculations of 
the induced interaction and noise in 3D bulk material. 
A comparison with the dipole-dipole spin interaction is 
given. 

Averaging the spin orbit Hamiltonian over the lo- 
calized electron's wave function, one obtains the spin- 
orbit coupling in the presence of magnetic field H — 



{Hx, Hy, Hz) in the form 



m,l — x^y,z 



(6.1) 



where /^b is the Bohr magneton. Here the tensor gmi 
is sensitive to lattice deformations. It was shown^^ that 
for the donor state which has tetrahedral symmetry, the 
Hamiltonian (|6.1[) yields the spin-deformation interaction 
of the form 

H = A/iB [SxxO-xHx + EyyayHy + e^^zCTzHz + A((T • H)/3] 
-I- B^B [sxy WxHy + UyHj,) -I- c.p.] . (6.2) 

Here c.p. denotes cyclic permutations and A is the effec- 
tive dilatation. The tensor already includes averaging 
of the strain with the gradient of the potential over the 
donor ground state wave function. 
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FIG. 8: The magnitude of the induced spin-spin interaction 
for a 3D Si-Ge type structure: The dominant interaction am- 
plitude, which is the same for the x and y spin components, is 
shown. The arrow indicates increasing b values for the curves 
shown, with b = 0.03, 0.09, 0.15, 0.21, 0.27, 0.33. The inset 
estimates the level of the noise (for 2kBTaB / Cs = 0.01): The 
bottom curve is ?7™(0), with m — x,y. For b < 0.4, the am- 
plitude 7?r'(d) can be comparable, and its values, calculated 
numerically for < r < 50, are shown as long as they exceed 
?7™(0), with the top curve corresponding to the maximum 
value, at r = 0. 



As before, let us assume that the separation d, as well 
as the magnetic field, are directed along the z-axis, for 
an illustrative calculation. Then the spin-deformation 
interaction Hamiltonian simplifies to 

H = A^J.b£zzO-zHz + B^B {SyzCTyHz + SzxCTxHz) ■ (6.3) 

In terms of the quantized phonon field, we havei*^ 



£-.,=^/(k)^ 

(6.4) 

where in the spherical donor ground state 
approximatio] 



.24.41 



/(k) 



(1 



(6.5) 



Here ob is half the effective Bohr radius of the donor 
ground state wave function. In an actual Si or Ge crys- 
tal, donor states are more complicated and include cor- 
rections due to the symmetry of the crystal matrix in- 
cluding the fast Bloch-function oscillations. However, the 
wave function of the donor electrons in our case is spread 
over several atomic dimensions (see below). Therefore, 
it suffices to consider "envelope" quantities. Thus, the 
spin-phonon Hamiltonian (|1.21 II. 3p coupling constants 
will be taken in the form 



{^k,zkm + Ck.mfcz) , (6.6) 



where — Dy — B/iB-ffz and = Afi^Hz- 
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Let us first consider a ID cliannel geometry along the 
z direction. This wiU give an example of an Ohmic bath 
model discussed towards the end of Section HIl In ID 
channel the boundaries'*^ can approximately quantize the 
spectrum of phonons along x and y, depleting the density 
of states except at certain resonant values. Therefore, the 
low-frequency effects, including the induced coupling and 
quantum noise, will become effectively onedimensional, 
especially if the effective gap due to the confinement is of 
the order of As mentioned earlier, this frequency cut- 
off comes from Ea. (l6.5p . namely, it is due to the bound 
electron wave function localization. A channel of width 
comparable to ~ ob will be required. This, however, may 
be difficult to achieve in bulk Si or Ge with the present- 
day technology. Other systems may offer more immedi- 
ately available ID geometries for testing similar theories, 
for instance, carbon nanotubes, chains of ionized atoms 
suspended in ion trapSj^^*^"^ etc. In our case, the longi- 
tudinal acoustic (LA, ||) phonons of the one-dimensional 
spectrum will account for the j| component of the cou- 
pling, whereas the transverse acoustic (TA, _L) phonons 
will affect only the x and y spin projections. 

One can show that the contributions of the crossprod- 
ucts of coupling constants, g]^^{g^^)* with m 7^ m', to 
quantities of interest vanish. The diagonal combinations 
are 



\9kJ 



\9kf 



AV|i?.' 



(6.7) 



\9l 



With our usual assumption for the low-frequency disper- 
sion relations ~ c\\kz and LOk^.± ~ c±kz, the ex- 
pressions (|6.7p lead to the Ohmic bath model discussed 
at the end of Section |TT1 The shape of the frequency 
cutoff resulting from Ea. (j6.7p is not exponential. How- 
ever, to estimate the magnitude of the interaction and 
noise one can utilize the results obtained earlier for the 
ID Ohmic bath model with exponential cutoff. The cou- 
pling constants in Eqs. (|2.17p and (|2.16p should, then, be 
taken as 



(6.8) 



X y ^ PB-"z 

= = 8^' 

where S* is a cross section of the channel, and the cutoff 
is ujc — > C||/aB for the z component, and lUc — > c^/ob 
for the X and y components. Considering Si as an exam- 
ple, we arrive at an approximately adiabatic Hamiltonian 
(ckf 3> ai'^) with Ohmic-type coupling. The dynamics of 
the concurrence, then, is qualitatively similar to the one 
shown on Fig. 2, with the peak entanglement ^ 0.4 as 
well. The coupling constant af , however, is significantly 
smaller due to the weakness of the spin-orbit coupling 



of P-impurity electrons in Si, which results in low mag- 
nitude of the induced interaction (and the noise due to 
the same environment), and slightly modifies the shape 
of the concurrence. 

In the 3D geometry, let us consider for simplicity only 
the LA phonon branch, ^ k/|k|, and assume an 
isotropic dispersion uj^,^ = |k|. The expression for the 
coupling constants is then 



Dr. 



(1 



(6.9) 



The cross terms, with m 7^ n, of the correlation func- 
tions TrB [Xi^X'^{t)pB], see Eq. (|X3l) in Appendix El 
depend on the combination 9^^ig^ ^)* , which is always 
an odd function of one of the projections of the wave vec- 
tor. The nondiagonal terms thus vanish, as mentioned in 
Appendix |^ 

Integrating Eas. ((2?8)) and ((2J)) with Eq. ((6J)) . see Ap- 
pendix [BJ one can demonstrate that dccoherence is dom- 
inated by the individual noise terms for each spin, with 
the typical amplitude 



15 (1 + 62)^ 



■ coth ■ 



2kBT' 



(6.10) 



where b = Aa^/cs and Cb = B'^p%H'^/ [16tt^ pha^cl) . 
The interaction amplitude x™(d) and, therefore, the in- 
duced spin-spin interaction, has inverse-square power-law 
asymptotic form for the x and y spin components, with a 
superimposed oscillation, and inverse-fifth-power-law de- 
cay for the z spin components 



Hint 



J2 2xTid)al^al 

^ — 

2b sin br 



(6.11) 



(1 + 62)4 ^2 



1.^2 



3847r^CA^cr> 



Here r = |d|/aB and Ca — A^Cb/B^. At small distances 
the interaction is regular and the amplitudes converge to 
constant values, see Figure [S] The complete expressions 
for xj^(d) and ?7™(d) are given in Appendix IB] 

In Figure [51 we plot the amplitudes of the induced 
spin-spin interaction (16.111 IB.4[ IB.7P and noise for dif- 
ferent values of the spin-spin separation and b, for elec- 
tron impurity spins in 3D Si-Ge type structures. The 
value of b can be controlled via the applied magnetic 
field, b — PBH^g* o-b/ Cg. The temperature dependence of 
the noise is insignificant provided 2kBT j ^ <^ 1. 

A typical value'^-'^'^ of the effective Bohr radius 
in Si for the P-donor-electron ground state wave func- 
tion is 2aB = 2.0 nm. The crystal lattice density is 
p = 2.3 X lO^kg/m^, and the g-factor g* = 1.98. For 
an order-of-magnitude estimate, we take a typical value 
of the phonon group velocity, Cs = 0.93 x 10* m/s. The 



spin-orbit coupling constants in Si ar& 
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FIG. 9: The magnitudes, measured in units of s~^, of the 
induced spin-spin interaction, the EM couphng strength, and 
the level of noise for P-impurity electron spins in Ge. Here 
Hz = 3 X 10'' G, and low temperature, 2kBT/A <C 1, was 
assumed. 



and B'^ « 10 ^. The resulting interaction constants in 
Eas. (l6T0l[6lI]) are Cb = 7.8 s^^ and Ca = IO^Cb. In 
the Ge lattice, the spin-orbit coupling is dominated by 
the non-diagonal terms^i^^i^ « and « 10^. 
The other parameters are 2aB = 4.0 nm, p = 5.3 x 
lO^kg/m^, cs = 5.37 x lO^m/s, and g* = 1.56. This 
results in a much stronger transverse component interac- 
tion, Cb ~ 1.3 X 10^ s^^ and Ca ~ 0. In both cases the 
magnetic field was taken Hz = 3 x 10^ G. In the above 
estimations, one could use various experimentally sug- 
gested values for the parameters, such as, for instance, 
qb- This will not affect the results significantly. 

As mentioned is Section |TT1 the obtained interaction 
(|6.1ip is always accompanied by noise coming from the 
same source, as well as possibly by other, direct inter- 
actions of the spins. When the electron wave functions 
overlap is negligible, the dominant direct interaction will 
be the dipole-dipole one 



An 



IdP 



(6.12) 



The comparison of the two interactions and noise is 
shown in Figure[9l We plot the magnitude of the effective 
induced interaction (|6.11l lK4l IB.7|) . the electromagnetic 
interaction (|6.12p . measured by T^em = A^oA^B/^^ldp, 
and a measure of the level of noise, for P-donor electron 
spins in Ge. It transpires that the induced interaction can 
be considerable as compared to the electromagnetic spin- 
spin coupling. However the overall coherence-to-noise ra- 
tio is quite poor for Ge. In Si, the level of noise is lower 
as compared to the induced interaction. However, the 
overall amplitude of the induced terms compares less fa- 
vorably with the electromagnetic coupling. 

In conclusion, we have studied the induced indirect ex- 
change interaction due to a bosonic bath which also in- 
troduces quantum noise. We demonstrated that it can 
create substantial two-spin entanglement. For an ap- 
propriate choice of the system parameters, specifically. 



the spin-spin separation, for low enough temperatures 
this entanglement can be maintained and the system can 
evolve approximately coherently for many cycles of its in- 
ternal dynamics. For larger times, the quantum noise ef- 
fects will eventually dominate and the entanglement will 
be erased. 

Estimates for P-impurity electron spins in Si and Ge 
structures have demonstrated that the induced interac- 
tion can be comparable to the dipole-dipole spin inter- 
action. One can also infer that this phonon-mediated 
interaction in the bulk (3D) Ge is not very effective to 
be used for quantum computing purposes, i.e., to entan- 
gle qubits. This is due to poor coherence to noise ratio. 
Therefore, the use of Si may be favored, despite the fact 
that it has weaker spin-phonon coupling. Indeed, the 
noise amplitudes for Si are significantly smaller then the 
induced exchange interaction, where the latter is domi- 
nated by the adiabatic term. The situation is expected 
to be further improved for systems with reduced dimen- 
sionality for phonon propagation. 
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APPENDIX A: DERIVATION STEPS FOR THE 
INDUCED INTERACTION AND NOISE 

An important assumption required for the validity of 
the Markovian approach concerns the time scale of the 
decay of the bath correlations introduced in Eqs. (|A.2[) 
and (|A.3[) below. By constantly resetting the bath to 
thermal, one implies that this time scale is significantly 
shorter than the dynamical system time scales of interest. 
The Markovian treatment, considered here, is, therefore, 
valid at all but very short times. The short-time dynam- 
ics, for the time scales down to order l/o^c, requires a 
different approachj^ii^i^' Note that one usually assumes 
that cjc ^ A. 

Here we review some of the steps that lead from 
the Markovian equation for the density matrix (12. 2p . to 
the expressions for the induced interaction and quan- 
tum noise. Substituting Eqs. pT^ and in Ea. (l^ . 
one can represent the integrand, TrB[HsB,[HsB{t' — 
t)i PBPs{t)\\Ps{t), as a summation over i,j,m,n, of the 
following expression: 

Ttb {XiMf - t)PB) «it' - t)ps{t) 
-Ttb {XlpBXl,{t' - t)) aMtWn{t' - t) (A.l) 
-Ti'B [Xl^if - t)pBXl) ai^it' - t)ps{t)al 
+Ttb (pBXUf - t)Xl) ps{t)al{t' - t)ai,. 
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Here crj,(t) = e'-f^^Vj^e-^-f^^*. All the terms in Eq. lfOI) 
involve the correlation functions 

C,nn ((1 - t) = TrBXi^Xl^{t)pB. (A.2) 

where X^^{t) — e*^^*X^e^'^^*. The explicit expression 
for these functions can then be obtained from Eas. (|1.3p 
and (|A.2|) . 

oo 

C™„(d,i) = — 3 ^ / dk\g^^gZ^\ fisinu;k,ct (A.3) 
(27r) J L 

^ —OO 

+ coth COS wk.et] COS (k • d + (/-C^^ - cj)]:^^) , 

where (/){^ ^ is a possible phase of the coupling constants 
g^^, which is not present in most cases. The coupling 
constants are examined in detail in Sections |V] and IVIl 
and explicit model expressions are given. Presently, we 
only comment that in many cases the resulting matrix 
of the correlation functions Cmn(d,i) is diagonal, which 
simplifies calculations, as illustrated in Section IVll 

The summation of Ea. (|A.l[ ) over i,i,m,n is further 
simplified by noting that C„m(d,i) ~ C^„(d, — t), and 
writing cr^(i) explicitly as ct* for n = z, and cr^ cos Ai + 
\[c\,cyl^ sin Ai for n ^ z. For a diagonal Cm.n, we then 
get the amplitude expressions 



and 



u 

j rftTma„m(d, t') cos(At') = xr(d), 



(A.4) 



U 

j diTmC„™(d,<')sin(At') 



xr(d), (A.5) 



U 

j di'ReC™„(d,t')cos(At') = 77;"(d), 

— oo 


j dt'ReC„™(d,t')sin(At') = r;,"(d). 

— OO 

This finally leads to Eqs.lg^ and ((23|) . 



(A.6) 
(A.7) 



APPENDIX B: INTERACTION AND NOISE 
AMPLITUDES FOR Si-Ge TYPE SPIN-ORBIT 
COUPLING 



In Eqs. ([2^ and jM]), with Eq. ljOl) . there is a com- 
mon angular part 

r r 

Im{kd) = dip smed0^^cos{kdcose), (B.l) 



which gives 



Ix,y{kd) = 47rRc 



12 - Pd'' 



ikd 



47rRei 



Pd^ 
5Pd^ - 12 



(B.2) 



,ikd 



I^{kd) = 167rRe 
— 47r Re i 



Pd^ ^ ifcrf 

k^d^ 

Pd^ - 12Pd^ + 24 
k^ ' 



(B.3) 



ikd 



Here d — |d|. The remaining integral in Eq.i 

oo 

Dl r Ix,y{kd) CPdk 



Xc'nd) = 



2p{2^fcsJ (l + a2fc2)^c2/c2- A2' 



(B.4) 

can be evaluated along a contour in the upper complex 
plane. The integration contour includes two simple poles, 
at fc = ±A/cs, which have to be taken with weight 1/2 
due to principal value integration. Also included is the 
pole at fc = i/flB, of order four, and a simple pole at = 
(with weight 1/2). The latter pole is for the second term 
in Eq. (|R^ only. 

The pole at fc = i/a^ yields exponentially decaying 
terms ~ exp(— d/ae). At large rf, the asymptotic behav- 
ior is controlled by the poles at fc = ztA/c^, 



XVI ^\ A 2^ 26 sin6r 



(B.5) 



where h = l^aB/cs and r — |d|/aB. The complete ex- 
pression can be easily obtained from Eq. (|B.4[) . One can 
also note that at d ^ 0, Eq. (jB.4p is 



X^'^d) = -APCb^-^^^^^^^ + 0{t^). (B.6) 
240(1 -^52)^ 

Along the same contour, the z component of Eq. (|2.8p . 

oo 

Dl, f Ix,y{kd)Pdk 



X^(d) 



(B.7) 



has only two poles, at fc = [order 2 for the first term, 
and also of order 1 or 3 for the second term in /^(fcd)], 
and at fc = i/aB (of order 4). The pole at fc = is to be 
taken with weight 1/2 and gives the l/r^ asymptotic. 



X^(d) 



3847r^CA-r, 



while for 



one obtains 



x:(d) = -C^- + 0(r2) 



(B., 



(B.9) 



Substituting Eqs. (|B.5| ) and (jB.Sp in iJint, i-e., in the sec- 
ond term in Eq. (|2.4p . one obtains Eq. (|6.1ip . 

By using Eqs. (|B.ip and (|B.2p . for Eq. (l^ we obtain 



€''"{<^) = ^Cb ^ 



(1 + 52)' 



-Ix.yirb) coth 



A 



2fcBT' 



(B.IO) 



The decoherence processes are dominated by the local 
noise terms, e.g., r7|?'^(0). Noting that lx,y{0) = one 



12 



obtains (|6.10p . One also finds that ry^(d) = 0. For low 
temperatures, ^k^T ^ '^'^^ coth 2k^T ^ ^ ai,nd, 
therefore, xfW «^ vTW, see Eqs.^J^ and (l2l|). Note 
that the expressions involving xT i^) '^^^^ be neglected in 
(|2.4p , since the corrections they introduce to the induced 



interaction have the same magnitude as the noise ampli- 
tudes. The function f?™(d) is often comparable to »7™(d) 
for Eq. (|6.9p with the parameters used in Sect ion IVTl This 
amplitude is calculated numerically in Figure [S] 
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